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ABSTRACT 

We revisit a core prediction of the disc instability model (DIM) applied to X-ray bina- 
ries. The model predicts the existence of a critical mass transfer rate, which depends on 
disc size, separating transient and persistent systems. We therefore selected a sample 
of 52 persistent and transient neutron star and black hole X-ray binaries and verified 
if observed persistent (transient) systems do lie in the appropriate stable (unstable) 
region of parameter space predicted by the model. We find that, despite the significant 
uncertainties inherent to these kinds of studies, the data are in very good agreement 
with the theoretical expectations. We then discuss some individual cases that do not 
clearly fit into this main conclusion. Finally, we introduce the transientness parameter 
as a measure of the activity of a source and show a clear trend of the average outburst 
recurrence time to decrease with transientness in agreement with the DIM predictions. 
We therefore conclude that, despite difficulties in reproducing the complex details of 
the lightcurves, the DIM succeeds to explain the global behaviour of X-ray binaries 
averaged over a long enough period of time. 

Key words: accretion, accretion discs - black hole physics - instabilities - methods: 
observational - X-rays: binaries. 
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1 INTRODUCTION 

Low mass X-ray binaries (LMXB) are bright X-ray sources 
consisting of an accreting black hole/neutron star primary 
and a Roche lobe filling, low-mass secondary star. Matter is 
transferred from the secondary star via the inner-Lagrange 
(LI) point and forms an accretion disc around the primary. 
In some cases the accretion disc undergoes sporadic out- 
bursts which are thought to be triggered by a thermal- 
viscous instability, resulting in an increased mass accretion 
rate onto the primary and a rapid X- ray brightening. This 
instability is thought to arise when the disc effective tem- 
perature is close to 10 4 K, enough for hydrogen to become 
partially ioni zed and for opacities to depend strongly on 
temperature (|Mever fc Mever-Hofmeisterlll98ll h The model 
describing the global properties of this instability applied to 
accretion discs is referr ed to as the Disc Instability Model 
(DIM; See lLasota|[200ll . for a review of the model). 

The DIM has been shown to reproduce the broad ob- 
servational properties of black holes and neutron stars out- 
bursts when it is modified to include irradiation heating and 
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a transition to a radiatively ineffic ient accretion flow below a 
given mass transfer rate ( King fc Ritter|[T998l ; iDubus et al.l 
bOQlf ). However, several features of outburst lightcurves re- 
main unexplained (glitches, plateaux, rebrightenings, etc.) 
and the model rests on critical assumptions for angular mo- 
mentum transport (alpha viscosity) and unknown physics 
(efficient /inefficient disc transition). Hence, it is difficult to 
use outburst lightcurves to provide definitive tests of the 
DIM. In contrast, one of the core predictions of the DIM, 
which depends little on the details of angular momentum 
transport, is the existence of a critical mass transfer rate that 
depends on disc size, separating stable and unstable discs. A 
basic test is therefore to verify if observed persistent (tran- 
sient) systems do lie in the appropriate stable (unstable) 
region of parameter space predicted by the DIM - modulo 
the effects of irradiation. 



van Paradiisj (|l996l ) carried out this test for the first 
time and showed that the DIM could indeed reproduce the 
transient-persistent dichotomy in X-ray binaries (XRBs) as 
long as irradiation heating is taken into account. The X-ray 
emission coming from the inner regions can indeed heat the 
outer parts of the accretion disc above the hydrogen ioniza- 
tion temperature and therefore stabilise the flow. Evidence 
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for the influence of irradiation comes mainly from the higher 
optical to X-ray luminosity ratio observed in XRBs com- 
pared to cataclysmic variables (CVs) which should be less 
affected by i rradiation due to smaller discs a nd accretion 
luminosities (Ivan Paradiis fc McClintocldll994h . 

Despite the li mited numbe r of so urces available at that 
time, the work of Ivan Paradiisl l| 19961 ) represents one of the 
most fundamental test of the applicability of the DIM to 
LMXBs. Fifteen years later, the number of known XRBs 
and the amount of data available to study them have in- 
creased significantly. We thus believe that it is time to re- 
visit and refine this test, which is t he aim of this paper. We 
also point out the recent work bv lJaniuk fe Czernvl (|201ll ) 
which discusses theoretical aspects of the hydrogen ioniza- 
tion and radiation pressure instabilities and present obser- 
vational constraints using a sample of Galactic black hole 
binaries. 

This paper is organised as follows. In Sect. 2, we detail 
the main assumptions and equations leading to the relation 
between the critical mass accretion rate and the orbital pe- 
riod on which the test rely. Sect. 3 presents our sample of 
sources and describes the data analysis and the method we 
used to estimate the mass transfer rate from the companion 
star. We present the results in Sect. 4 and discuss the out- 
comes in Sect. 5. The last section summarises and concludes 
this work. 



2 A SIMPLE TEST OF THE DISC 
INSTABILITY MODEL 

An accretion disc is globally stable, and does not experience 
outbursts, whenever all its annuli are stable. The disc can be 
in a hot, stable equilibrium if the local accretion rate M(R) 
is higher than the critical accretion rate of the hot, ionized 
state, M cr it (R), for all radii R in the disc. Since M cr it(-R) 
increases with radius, this condition can be fulfilled if the 
mass transfer rate from the companion star, M ex t, is higher 
than the the critical accretion rate at the outer radius of the 
disc M cr it (Rout)- This is the standard scenario for persistent 
X-ray binaries. 

However if M ext < M CI it(R ou t), the disc enters a limit 
cycle behaviour, oscillating between the ionized state (out- 
burst) and the neutral state (quiescence), i.e. the system is 
transient. 

When the effect of irradiation is neglected, the critical 
mass accretion rate above which a ring of matter located 
at a rad ius R is in the stable hot sta te can be expressed as 
follows jLasota. Dubus fc Kruk|[2"008l . appendix A): 



where C = 10 C_3 is the irradiation constant defined by 
the irradiation flux: 



A> q r\*7 1n l5 — 0.01 D 2.64 ,,-0.89 -1 

M C Ht = 8.07 x 10 a 01 R 10 M 1 gs 



(1) 



where Mi is the mass of the compact object in solar units, 
a = 0.1 ao.i is the viscosity parameter and R — Rio 10 10 cm. 
When irradiat ion is taken into account, the c ritical accretion 
rate becomes l|Lasota. Dubus fc Krukll2008l . appendix A): 



T 4 r Mc 2 



(3) 



The above expression of M cr it has been derived by in- 
vestigating the local vertical structure of an irradiated disc 
where X-ray heating is limited to a thin layer above the disc 
and appears in th e boundary condition at the surface (see 
iDubus et al]|l999l . for details): 



ry-i4 /jt4 - rp4 

■*■ surf ■*■ eff T -* irr 



(4) 



where T BUT { is the surface temperature and T er r the effective 
temperature in the absence of irradiation. 

As mentioned previously, the overall disc will be sta- 
ble if the mass transfer rate from the donor star M ex t is 
higher than the critical accretion rate (Eq. [T] or at the 
outer radius Rout- For a given mass ratio q = M2/M1, one 
can compute the ratio R out /a between the outer radius and 
the orbital separ ation following the calculations given in e.g. 
iPaczvnskl (|l977l ). The orbital separation can then be written 
as a function of the orbital period P (in hr) using Kepler's 
law: 



a = 3.53 x 10 10 Mi /3 (1 + q) 1/3 f£ /3 cm 



(5) 



Consequently, for a given set of parameters a, C, Mi and q, 
we can derive a relation between the critical accretion rate 
at the outer radius of the disc and the orbital period. 

The test therefore consists in estimating the mass trans- 
fer rate of a sample of persistent and transient XRBs (black 
holes and neutron stars) with known orbital periods and 
compare their position in a M — P diagram with the tran- 
sient/persistent separation predicted by the M cr it(-P) rela- 
tion derived from the DIM. 

Note that the irradiated case (Eq. [2} contains a free 
parameter, the irradiation constant C. The latter is a mea- 
sure of the fraction of the X-ray luminosity intercepted and 
thermalised in the upper layers of the outer disc and, as 
such, contains information on the irradiation geometry, ra- 
diative efficiency, X-ray albedo, X-ray spectrum etc. It is 
thus a poorly constrained quantity. Its con stan cy can also 
be qu estioned. lEsin. Lasota fc Hvnes! (|2000T ) and lEsin et al.l 
(2000) present some evidence of changes of C during the out- 
bursts of GRO J1655-40, A0620-00 and GRS 1124-68. The 
rece nt results on accret ion disc winds in black hole XRBs 
from iPonti et all (|2012h also suggest a significant variation 
of the influence of irradiation during outbursts. Nonetheless, 
for simplicity, we restricted ourselves to constant C in this 
study and started by assuming a value of 10~ 3 found to be 
compatible with the observed optical magnit udes and sta- 
bility properties of persistent X-ray binaries (|Dubus et al.l 
1999). This value is also consistent with the results from, 
e.g., Ide Jong et ail l| 19961 ) or iGierliriski etafl (|2009l ) who 
found C ~ 10~ 2 — 10 -3 by comparing irradiated disc models 
to optical/UV and X-ray data. 



M C rit = 9.5 x 10 1 * CZ3 

X Mi -°-<*+0.081 OgC - 3gs -l 



0.36 0.04+0.01 log C_ 3 r.2.39-0.10 logC_ 3 
H 



10 



(2) 



3 DATA SELECTION AND ANALYSIS 

To test the validity of the stability criteria (Eq.[T]and[2}, we 
selected a sample of 52 XRBs with known orbital periods, 
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reasonably well-sampled light curves and for which we think 
that accretion proceeds via a disc. Twenty three of these 
contain a black hole (BH) or a black hole candidate (BHC) 
and the other 29 host a neutron star (NS). This is more than 
twice the number of sources used in the original work by 
Ivan Paradiisl l| 19961 ). We classified the sources as BH or NS 
according to the estimates of the mass of the compact object 
and the detection of nuclear bursts. If the BH nature of 
the compact object is only suggested by the global spectral 
and temporal properties of the source, we classified it as a 
BHC. However, the exact nature of the compact object has 
a little influence on the results of the test since the Merit (P) 
relation derived from the DIM is very weakly dependant on 
the primary mass (see Sect. [4]). 

All the sources in our sample are considered as low-mass 
XRBs except the three BH systems Cyg X-l, LMC X-l and 
LMC X-3 and the three NS systems SMC X-l, LMC X-4 and 
Cen X-3. These objects likely host a high-mass companion 
star producing winds that are later accreted by the compact 
object. However, the X-ray spectro-temporal properties of 
these sources indicate the presence of an accretion disc fed 
by these focused winds. In some cases (e.g. Cen X-3), it is 
very likely that the companion star fills its Roche lobe and 
therefore transfers mass through the LI point in addition to 
the winds. We therefore included them in our sample taking 
into account that the size of the disc derived from the orbital 
period might be overestimated. 

The BH and NS samples contain respectively 6 and 17 
sources that we classified as persistent systems. To do so, 
we assumed that the X-ray sk y is well covered (~ 80 — 90%) 
since 1970 (IChen et al.l ll997h. If the source was discovered 
in the early 1970s or before and is active since then, we 
classified it as persistent. By "active", we mean here that it 
never went below the detection threshold of the instruments 
that observed it throughout the years. 

Note that two sources do not easily fit into this simplis- 
tic bimodal classification. The BHCs and microquasars IE 
1740.7-2942 and GRS 1758-258 have been discovered respec- 
tively by the Einstein satellite in 1984 dHertz fe Grindlavl 
1 1984 ) and the Granat satellite in 1990 (|Mandroul Il990h 
Both sources are commonly considered as persistent in the 
literature since they remained active since their discov- 
ery. No X-ray source was identified at the location of IE 
1740.7-2942 during previous observations of the Galactic 
cent re region by the sate llites Uhuru, Ariel V and SAS 3 
fsee lPavlinskv et al.lll992l and references therein). However, 
given the variability of the source and the sensitivity and 
angular resolution of these instruments, we cannot exclude 
that IE 1740.7-2942 was already active but not detected. 

In the case of GRS 1758-258, due to its proximity (~ 40 
arcmin) to the persistent neutron star Z-source GX 5-1, it is 
very likely that both sources were considered th e same before 
Gran at was able to identify them separately (|Levine et al.l 
1 19841 already suspected that the hard X-ray flux assigned 
previously to GX 5-1 may originate from another nearby X- 
ray source). Given that an X -ray source is already present 
in the second Uhuru catalog l|Giacconi et al.|[l972T ) at a lo- 
cation consistent with GRS 1758-258, it is possible that the 
source was already active before its identification. Bearing 
all these considerations in mind, at first instance, we will 
treat these sources as persistent 

We also stress that among the 6 persistent BH systems, 



3 of them are the high-mass XRBs mentioned previously and 
2 of them are the black hole candidate IE 1740.7-2942 and 
GRS 1758-258 for which, as we just warned, the persistent 
nature is unclear. The last one, 4U 1957+115, agrees well 
with our persistency criterion. The sourc e remained per- 
sisten tly in soft states since its discovery (|Wiinands et al.l 
2002) and recently provided very deep radio upper limits 
supporting the jet quenching p aradigm in soft states black 
hole XRBs l|Russell et al.ll201ll ). The X -ray properties of the 
source s trongly suggest th e system is hosting a black hole 
(see e.g. iNowak et al.|[2008T ) but there is no firm evidence on 
the nature of the compact object. The persistent BH sample 
should thus be considered with caution. 

Note that the proportion of sources with known orbital 
periods corresponds to ~ 50% of the known sources in the 
case of BH and about 30% in the case of NS systems. These 
proportions must be kept in mind with respect to the results 
of this study. The sample we have selected might also be 
biased toward high inclination sources for which it is easier 
to determine the orbital period. 

For each source we selected, we list in Table[T]and[2]the 
mass, distance and orbital period we used for this study. 



3.1 Estimating the mass transfer rate 

The bulk of the work presented here was to estimate, for 
each source, the mass transfer rate from the donor star. For 
persistent sources, this task is fairly straightforward. If we 
assume that the accretion rate in the disc, averaged over 
time, is equal to M cxt , we can estimate the latter by mea- 
suring the average X-ray luminosity: 



L x = J7M ox tC 



(6) 



where c is the speed of light and n is the radiative efficiency 
of accretion. 

For transient systems, the standard procedure to es- 
timate Mext is based on the assumption that the systems 
mainly store mass during quiescence and accrete it suddenly 
during an outburst. The rate at which the mass accumulates 
in the disc, M accU m, can be estimated from the integrated 
X-ray luminosity from the system during an outburst, AE, 
and from the outburst recurrence time, At: 



AE 
Atnc 2 



(7) 



If the mass accretion rate during quiescence is negligi- 
ble compared to the mass transfer rate from the sec- 
ondary, we can then use M accum as a proxy for M ext . 
This assumption ignores the possibility of a highly ineffi- 
cient accretion flow which could hide a substantial accre- 
tion of mass in quiescence jNarayan. Barret fc McClintockl 
1 19971 : iMenou. Naravan fe Lasotal Il999i ). Furthermore, esti- 
mating Af acC um from X-ray emission alone does not take 
into account mass loss thr ough winds and jets (e.g., 
iFender. Gallo fc Jonkerl 120031 1. Consequently, for both rea- 
sons, our estimates of M ex t should be considered as lower 
limits to the actual mass transfer rate. 

For transient BH systems, the radiative efficiency may 
vary during an outburst. It is generally thought that the 
accretion flow is radiatively inefficient during faint hard state 
with r\ oc M while it is radiatively efficient during soft states 
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Table 1. Black Holes X-ray Binaries Parameters 



Source 


Type 


Af(M ) 


D(kpc) 


Porb(hr) 


Notes 


GRO J0422+32 


T 


3.97+0.95 


2.75 ±0.25 


5.092 




A 0620-00 


T 


11.0+1.9 


1.2 ± 0.4 


7.75 




GRS 1009-45 


T 


> 3.9 


5.7 ± 0.7 


6.84 




XTE J1118+480 


T 


8.53 


1.8 ± 0.6 


4.08 




GS 1124-684 


T 


4-11 


5.5 ± 1.0 


10.38 




GS 1354-64 


T 


> 5.75 


^25 


61.07 


(1) 


4U 1543-47 


T 


8.4-10.4 


7.5 + 0.5 


26.8 




XTE J1550-564 


T 


10.5+1.0 


5.3 + 2.3 


37.25 




XTE J1650-500 


T 


> 2.73 


2.6 ± 0.7 


7.63 


(2) 


GRO J1655-40 


T 


6.3+0.5 


3.2 + 0.2 


62.88 


(3) 


MAXI J 1659-352 


T, BHC 


3.6-8.0 


8.6 


2.4 


(4) 


GX 339-4 


T 


> 5.3 


~8 


42.14 


(5) 


4U 1705-250 


T 


> 4.73 


8.6 + 2.0 


12.54 




Swift J1753.5-0127 


T, BHC 




~8 


3.2 


(6) 


GRS 1915+105 


T 


14+4 


9.0 + 3.0 


739.2 


(7) 


GS 2000+25 


T 


4.8-14.4 


2.7 + 0.7 


8.26 




V404 Cyg (GS 2023+338) 


T 


12+2 


2.39 + 0.14 


155.4 


(8) 


IE 1740.7-2942 


P?, BHC 




~8.5 


305.5 


(9) 


GRS 1758-258 


P?, BHC 




~8.5 


442.8 


(9) 


4U 1957+115 


P, BHC 




3= 7 


9.33 


(10) 


Cyg X-l 


P, HMXB 


14.8 ± 1.0 


1.86 + 0.12 


134.4 


(11) 


LMC X-l 


P, HMXB 


10.91 ± 1.41 


48.10 + 2.22 


93.6 


(12) 


LMC X-3 


P, HMXB 


11.1 ± 1.4 


52.0 + 0.6 


40.8 


(13) 



Col. (Source): the source name. (Type): "T" and "P" indicate transient and persistent systems respectively. "HMXB" indicates a high 
mass X-ray binaries. (M): the mass of the co mpact object. (D): the distance to the source. (P or b): t he or bital period of the system. 
Unless otherwise stated, parameters are fromlLi u , van Paradiis fe van den Heuvel J2005I. l2006t 120071) an d I Jonker fc Nelemansl ll2004h 
NOTE S: (1) Lower limit on the distance from ICasares^T*aT j2009h . (2) Dista nce from iHoman et al l ll2006h. (3) Note that iFoellmil 



i20091 ) argues for an upper limit to the distance of 2 kpc. (4) Parameters fro m iKuulkers^^^^""!^^^) and lYamaoka et al] l l2012h (5) 
Distance from IZdziarski et al.l [|2004h . (6) Par ameters from IZurita et al l J2008h and ICadolle Bel et al] (|2007f ). (7) O rbital period from 
iNeil, Bailvn fc C obb (2007) and distance from I Fender et al I jl99Sh and lChapuis fc Corbell i2004h . (8) Distance from iMiller- Jones et al.l 
l|2009jV_ (9) No estimat e of the mass of th e compact object . We adopted a Galacti c center distan ce for the so urce. Orbital period 
fromlSmith et al.l j2002lV (10) Distance fromlYao et al] j200St) andlNowak et all d200Sf) (see al solNowak et al.ll2012T) . orbital period from 
Thorstensenl il987ft andlBavless et al.l l|201ll ). ( 11) Mass from Orosz et al.l l|201ll l. dis tance fromlReid et aljfcOlll ) and orbital period from 
Brock sopp et al] il999l ). (12) Parameters from Orosz et al.l j2009l'l . (13) Mass from Gicrliriski et al. I ll2001r i, distance from Idi Benedetto] 



and orbital period from Hutchi ngs et al. I l|2003r i 



with 77 ~ 0.1. In addition, the total X-ray luminosity usually 
remains relatively constant (within a factor of ~ 2, e.g., 
IZhang et al]ll997ft during the hard to soft state transition 
which does not suggest a drastic change in the radiative 
efficiency. Therefore, during the bright hard state phases, 77 
should not be significantly lower than 0.1. Consequently, we 
adopted the following prescription for 77 when deriving M ex t 
for a transient BH XRB: 

v = 0.1 ( ™ ) ; for L x < l%L Edd 

V0.01M E dd/ (8) 

T) = 0.1 ; for L x > l%L E dd 

where the Eddington accretion rate is denned by L E dd = 
0.1Af E ddC . Accretion onto neutron stars is assumed to be 
radiatively efficient due to the presence of a solid surface on 
which the accreted material can eventually radiate its re- 
maining gravitational energy. We thus used 77 = const = 0.1 
to derive M ex t for NS systems. Note that we have chosen 
0.1 as the maximum value for 77, although the radiative ef- 
ficiency can reach higher values for matter rotating around 
highly spinning black holes. However, 77 departs significantly 
from 0.1 on ly for very high spin values (e.g., 77 ~ 0.2 for a 
spin of 0.9, iBardeen et al.|[l972i ), therefore, this shouldn't 
affect our results significantly for most of our sources. 



To obtain reliable measures of AE, we used all ob- 
servations of the binaries we selected which were publicly 
available in the archive of the Rossi X-ray Timing Explorer 
(RXTE) satellite. The data reduction and model fitting 
scheme is detailed in the next section. Some objects we se- 
lected have been inactive since the launch of RXTE in De- 
cember 1995. We thus used estimates of either AE or M ex t 
found in the literature that we corrected for updated dis- 
tance if necessary. From model fitting to the X-ray spectra, 
we produced X-ray light-curves in the 0.1-200 keV energy 
band for each sources. Then, we fit the light-curves with 
combinations of Gaussians or polynomial functions to inte- 
grate the luminosity and derive the energy liberated during 
each outburst. In the case where an outburst was covered 
by the All Sky Monitor (ASM) but only partially covered 
by the Proportional Counter Array (PCA) pointed observa- 
tions (typically, the rising phase was missing), we extrapo- 
lated the function used to fit the PCA lightcurve in order 
to match the shape of the ASM lightcurve. If the outburst 
was entirely missed by the PCA, we estimated the total en- 
ergy liberated using the ASM lightcurve and the conversion 
of 1 ASM coun t /s to 7.7 x 10~ 10 ergs _1 cm~ 2 proposed by 
lin't Zand et~ai1 (|2007l ). 
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Table 2. Neutron Star X-ray Binaries Parameters 



Source 




Type 


Af(M ) 


D(kpc) 


Porb(hr) 


Notes 


Tf^R 700901 

ivjxx juuzyi-t-oyo'y: 




T 




9 ^ f\ 


9 A(\ 




EXO 0748-676 




T 




5.9-7.7 


3.82 




Cen X-4 




T 


1.3+0.6 


1.2 + 0.2 


15.10 




4U 1608-52 




T 




2.8-3.8 


12.89 




XTE J1710-281 




T 




12-16 


3.28 




1A 1744-361 




T 




<9 


1.62 




GRS 1747-312 




T 




9.5 ± 3.0 


12.36 




XTE J1751-305 




T 




~8.5 


0.71 




SAX J1808.4-3658 




T 


<2.27 


3.4-3.6 


2.01 




XTE J1814-338 




T 




8.0+1.6 


4.27 




Aql X-l 




T 




4.4-5.9 


18.95 




XTE J2123-058 




T 


1.5+0.3 


8.5+2.5 


5.96 




LMC X-2 




P 




50 + 2 


8.16 




4U 0614+091 




P 




2-3 


0.8 


(1) 


All 1^9^ f^9 




1 5 
1 




in 
~ 1U 


9 Q 




Sco X-l 




P 




2.8+0.3 


18.9 




4U 1636-536 




P 




3.7-4.9 


3.8 




GX 349+2 




P 




5.0 + 1.5 


22.5 


(2) 


Her X-l 




P 


1.1 ± 0.4 


6.6 ±0.4 


40.8 


(3) 


4U 1735-444 




P 


> 0.1 


8.0-10.8 


4.6 




4U 1746-370 




P 




11.0 


5.16 




4U 1820-303 




P 




7.6+0.4 


0.19 


(4) 


4U 1850-087 




P 




8.2+0.6 


0.34 


(4) 


4U 1916-053 




P 




8.9+1.3 


0.83 




4U 2129+12 (AC 211) 




P 




10.3 + 0.4 


5.96 


(4) 


Cyg X-2 




P 


1.78 ± 0.23 


7.2 + 1.1 


236.2 




SMC X-l 


P. 


HMXB 


1.04 ± 0.08 


60 


93.36 


(3) 


LMC X-4 


P. 


HMXB 


1.29 ± 0.05 


50 


33.79 


(3) 


Cen X-3 


P, 


HMXB 


1.49 ± 0.08 


9 ± 1 


50.16 


(3) 



Col. (Source): the source name. (Type): "T" and "P" indicate transient and persistent systems respectively. "HMXB" indicates a high 
mass X-ray binaries. (M): the mass of the co mpact object. (D): the distance to the source. (P or h): t he or bital period of the system. 
Unless otherwise stated, param eters are from [Liu, van Paradiis fc van den Heuve I J20051. 120061. l2007h and lJonker fe Nelemansf J2q04h. 
NOTES: (1) Orbit al period fromlin't Zand et al.1 l|2007l V (2) Distance from lChristian fe Swankl j 19971). (3) Mass from lRawls et al.l ^201lft . 
(4) Distance from lKuulkers et al.l (|2003| V 



Some objects underwent several outbursts which allows 
for a precise estimate of the recurrence time. Some other 
entered only once into an active phase. In that case we used 
a lower limit to At (leading to an upper limit to M ex t) by 
considering that t he X-ray sky is well covered since 1970 
l|Chen et al.lll997l ). Therefore, At > max(Ati 97 o, At2oi2), 
where Ati97o and Ai20i2 are the times elapsed from 1970 
to the outburst start date and from the outburst to 2012 
respectively. 

3.2 RXTE data reduction 

We reduced RXTE data using the HEAS0FT software package 
version 6.11 following the standard procedure described in 
the RXTE cookboolfl 

We extr acted spectra from the Proportional Counter 
Array (PCA; |jahoda et al.ll200rj) for which we only used the 
Proportional Counter Unit (PCU) 2. This is the only oper- 
ational unit across all the mission and is the best-calibrated 
detector out of the five PCUs. A systematic error of 1 per 
cent was added to all channels. We obtained hard X-ray 
spectra from the High Energy Timing Experiment (HEXTE) 

1 http:/ /heasarc. gsfc.nasa.gov/docs/xte/data_analysis. html 



using both cluster A and B until December 2005. After this 
date, due to problems in the rocking motion of Cluster A, 
we extracted spectra from Cluster B only. On the 14th of 
December 2009, Cluster B stopped rocking as well. From 
this date we thus used only PCA data in our analysis. 



3.3 Model fitting 

We performed simultaneous fit s to the PCA and HEXTE 
spectra in Xspec version 12.7.0 l|Arnaudlll996l ) using a float- 
ing normalisation constant to allow for cross-calibration un- 
certainties. Channels out of the 3 — 50 keV range for PCA 
spectra and out of the 20 — 200 keV range for HEXTE spec- 
tra were ignored from the fit. 

The main objective of the X-ray spectral analysis was to 
obtain a reasonable estimate of the bolometric luminosity. 
We have chosen the 0.1-200 keV band as a proxy for the 
bolometric luminosity. Note that the RXTE data do not 
constrain the models outside of the 3-200 keV range. The 
0.1-200 keV flux is thus obtained by extrapolating the best- 
fitting models using an energy response matrix extended 
down to 0.1 keV (using the extend command in Xspec). 

By choosing the 0.1-200 keV flux as a proxy for the 
bolometric flux, we are neglecting a fraction of the accretion 
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luminosity. To estimate the importance of this missing frac- 
tion, we simulated with Xspec an irradiated accretion disc 
with a comptonisation comp onent using the diskir model 
iGierliriski et alj|2008l . l2009h . Very briefly, this model is an 
extension of the diskbb model that includes comptonisation 
of the blackbody di sc's photons by a corona of hot electron s 
(based on NTHCOMp lzdziarski et al.lll996l ; IZvcki et aljl998h . 
as well as disc irradiation of both the inner and outer re- 
gions of the disc. We first simulated typical soft states and 
h ard states spectra of BH XRBs using the parameters given 
in lGierlinski et all |2009T l. We then calculated the ratio be- 
tween the 0.1-200 keV flux and the 10" 6 - 500 keV flux 
(chosen to represent the bolometric flux). We found that, in 
the less favourable case, the 0.1-200 keV flux represents 75 
% of the bolometric flux. 

However, since we fit the data above 3 keV only, we also 
need to consider the case where the disc was missed because 
its contribution above 3 keV was not significant enough (typ- 
ically when fcTdisc < 0.4 keV). Extrapolating our best fit- 
ting model down to 0.1 keV implies that we only considered 
the contribution of the comptonisation component in the 
0.1-3 keV band and therefore probably missed a significant 
fraction of the bolometric flux. To estimate the fraction we 
missed in these particular cases, we simulated a hard state 
spectrum using diskir with fcTdi sc = 0.4 keV and measured 
the bolometric flux. We then compared the latter with the 
0.1-200 keV flux obtained using a simple (powerlaw + high 
energy cutoff) model with the same normalisation and spec- 
tral parameters than the comptonisation component of the 
diskir spectrum. We found that the 0.1-200 keV in that 
case represents ~ 50% of the bolometric flux. Consequently, 
during observations where the disc's inner temperature is 
low, we might underestimate the bolometric flux by a fac- 
tor of 2. However, these situations occur essentially during 
the low luminosity states. For the purpose of estimating the 
mass transfer rate, we are mainly interested in the high lumi- 
nosity states where most of the mass is accreted. Therefore, 
our results shouldn't be affected significantly. 

Significant variation of the hydrogen column density 
may occur during outburst of transient XRB s (see e.g. 
ICabanac et al . 2009, but also iMiller et al]l2009h . However, 
the PCA response falls below 2 keV and doesn't allow the 
interstellar absorption to be constrained properly. We thus 
fixed the Nu to the commonly accepted values for the bi- 
naries. We used the p hotoelectric absorption model p habs 
using abundanc es of [ Wilms. Allen fc McCravl d2000h an d 
cross-sections of iBalucinska-Church fe McCammonl I l992h . 



Table 3. Black Holes X-ray Binaries 



Given the high number of spectra to model, we auto- 
mated the fitting proce dure. For BH sou rces, we followed 
the method described in iDunn et all l|201(f) for which we re- 
fer the reader for full details. F or NS systems, we followed 
iLin. Remillard fc Homanl (|2007T ) and adopted their hybrid 
model (blackbody + broken powerlaw for hard states and 
disc blackbody + blackbody + constrained broken power- 
law for soft states) to fit the spectra. A full description of the 
model and its performance again st several desirability cri- 
teria (e.g. Lx oc T 4 ) are given in lLin. Remillard fc Homanl 
l|2007h . 

Table [3] and [4] present our estimates of the mass transfer 
rate and the average outburst recurrence time (if relevant) 
for our BH and NS samples respectively. 



Source 


Afcxt 


At av 


M oxt /M cri t 




(gs- 1 ) 


(yr) 




GRO J0422+32 


< 1.7 x 10 15 (1) 


> 22 


< 2.6 x 10~ 2 


A0620-00 


3.4 x 10 15 (1) 


58 


2.7 x 10" 2 


GRS 1009-45 


< 1.7 x 10 16 (1) 


> 23 


< 0.16 


XTE J1118+480 


1.3 x 10 16 


5 


0.28 


GS 1124-684 


< 2.1 x 10 16 (1) 


> 21 


< 0.10 


GS 1354-64 


> 1.0 x 10 17 


10 


> 2.9 x 10" 2 


4U 1543-47 


7.1 x 10 ie 


10 


7.7 x 10" 2 


XTE J1550-564 


1.0 x 10 17 


20 


6.5 x 10~ 2 


XTE J1650-500 


< 6.9 x 10 15 


> 31 


< 5.5 x 10~ 2 


GRO J1655-40 


5.5 x 10 16 


10 


1.6 x 10" 2 


MAXI J1659-352 


8.9 x 10 15 


> 40 


< 0.44 


GX 339-4 


8.0 x 10 17 


2 


0.43 


4U 1705-250 


< 1.8 x 10 16 (1) 


> 24 


< 6.5 x 10~ 2 


Swift J1753. 5-0127 


< 7.2 x 10 16 


> 35 


< 2.3 


GRS 1915+105 


< 9.8 x 10 18 


> 22 


< 5.5 x 10~ 2 


GS 2000+25 


< 4.1 x 10 15 (1) 


> 23 


< 2.9 x 10~ 2 


V404 Cyg 


5.9 x 10 16 (1) 


25(2) 


3.9 x 10" 3 


IE 1740.7-2942 


3.4 x 10 18 




7.7 x 10~ 2 


GRS 1758-258 


4.2 x 10 18 




5.3 x 10~ 2 


4U 1957+115 


> 9.2 x 10 16 




> 0.54 


Cyg X-l 


8.9 x 10 17 




< 7.5 x 10~ 2 


LMC X-l 


4.1 x 10 18 




< 0.61 


LMC X-3 


4.2 x 10 18 




< 2.34 



(1) Estimated from 
(2) Recurrence time 



outburst fluences given in IChen et al.l 419971 1 
estimated from lRichterl il989T ) 



Table 4. Neutron Stars X-ray Binaries 



Source 


Moxt 


Atav 


M oxt /M crit 




(gs- 1 ) 


(y) 




IGR J00291+5934 


3.2 x 10 14 


3 


2.3 x 10~ 2 


EXO 0748-676 


< 2.8 x 10 16 


> 28 


< 1.0 


Cen X-4 


2.4 x 10 15 (1) 


10 


1.0 x 10~ 2 


4U 1608-52 


6.1 x 10 16 


0.56 


0.33 


XTE J1710-281 


< 1.0 x 10 16 


> 28 


< 0.86 


1A 1744-361 


< 1.1 x 10 16 


2 (2) 


< 1.5 


GRS 1747-312 


7.1 x 10 16 


0.37 


0.40 


XTE J1751-305 


3.8 x 10 14 


3.8 (3) 


0.20 


SAX J1808.4-3658 


1.1 x 10 15 


2 


0.11 


XTE J1814-338 


< 3.8 x 10 14 


> 33 


< 1.2 x 10~ 2 


Aql X-l 


3.8 x 10 16 


0.55 


0.11 


XTE J2123-058 


< 4.3 x 10 14 


> 28 


< 7.8 x 10~ 3 


LMC X-2 


2.4 x 10 18 




27 


4U 0614+091 


2.9 x 10 16 




13 


4U 1323-62 


6.1 x 10 16 




3.5 


Sco X-l 


1.8 x 10 18 




5.2 


4U 1636-536 


7.9 x 10 16 




2.9 


GX 349+2 


1.5 x 10 18 




3.4 


Her X-l 


8.0 x 10 17 




0.68 


4U 1735-444 


4.0 x 10 17 




11 


4U 1746-370 


2.6 x 10 17 




6.0 


4U 1820-303 


7.9 x 10 1B 




3.4 x 10 2 


4U 1850-087 


2.4 x 10 16 




41 


4U 1916-053 


8.2 x 10 16 




34 


4U 2129+12 


2.5 x 10 17 




1.6 


Cyg X-2 


1.9 x 10 18 




0.1 


SMC X-l 


5.5 x 10 18 




1.3 


LMC X-4 


6.7 x 10 18 




7.7 


Cen X-3 


9.4 x 10 17 




0.58 



(1) Estimated from outburst fluences given in lChen et al.l il997T ) 
(2) Recurrence time and accretion rat e estimated from the act iv- 
ity of the source since 2003. (3) From lMarkwardt et all J2002f> . 



3.4 Dominant sources of uncertainty in the 
estimate of M ex t 

The value of M ex t we derive using the method described 
above will of course depend on the estimate of the distance 
to the source. If we assume an error of 25% on the distance, 
this will introduce an uncertainty of a factor of 2 on the 
estimate of M ex t- However, one can calculate using Tables [1] 
and [5] that the average error on the distance for the sources 
of our sample is ~ 20% for the black holes and ~ 15% for 
the neutron stars. Therefore, an uncertainty of a factor of 2 
on the value of M GXt should be representative. 

Owing to the smaller projected area of the disk, edge- 
on systems appear less luminous than face-one systems for 
the same intrinsic luminosity. The observed apparent X-ray 
luminosity is reduced by a factor of ~ cosi, where i is the 
angle between the line of sight and the rotation axis of the 
disc. Inclination could also affect the observed X-ray lumi- 
nosity if e.g. the disc starts to become geometrically thick 
and self-screens part of its inner regions, and/or if an equa- 
torial wind from the disc becomes optically thick, so that 
the luminosity seen at high inclinations is lower than the 
intrinsic flux seen at low inclinations. If these effects are sig- 
nificant, the mass transfer rate we infer from the total X-ray 
luminosity would be underestimated. 

Errors introduced by the model fitting procedure are 
not easy to estimate given that we extrapolate the models 
out of the energy range covered by the data. After running 
several tests comparing fluxes obtained using extrapolated 
models from RXTE data only and fluxes obtained with mod- 
els constrained at low energies by Swift data, we estimated 
that a typical error of a factor of 2 on the outburst fluences 
(flux integrated over time) introduced by the fitting proce- 
dure should be reasonable. 

An additional uncertainty is introduced by the efficiency 
factor T) with which gravitational energy is transformed to 
radiation and how representative the flux is for the mass 
accretion rate. However, the main contribution to the total 
energy liberated during an outburst arises mostly from the 
brightest phases where we assume rj — 0.1. This value is 
usually considered as an upper limit on the actual radiative 
efficiency for most of the sources, which means a lower limit 
on Mext- Therefore, this strengthen our previous caveat on 
our estimates of the mass transfer rate. 

Estimating the mass transfer rate in transient systems 
partly relies on how precisely we are able to estimate the 
recurrence time, which in turn depends on how we define 
quiescence and outburst. As a basic definition, we consid- 
ered that an outburst started (ended) when the source went 
above (below) the detection level of the ASM or the PCA 
depending on which one was the first (last) to detect it. 
Although very crude and instrument-dependent, this is the 
best definition we can use since regular monitoring of sources 
at very low luminosity is not yet possible. Since this defini- 
tion of quiescence will systematically overestimate the inter- 
outburst time, this will lead to a lower limit on the actual 
mass transfer rate, strengthening again our previous caveat. 
For sources with decade-long quiescent period, an overesti- 
mate of the order of weeks or even month on the recurrence 
time will have a negligible impact on the estimate of M ex t- 
For sources with short recurrence time (At < 2 yrs, e.g., GX 
339-4 for the black holes or Aql X-l and 4U 1608-52 for the 
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Table 5. Impact of the main sources of uncertainty 



Source of uncertainty 


Impact on M ex t 


High quiescent M 


Underestimated 


Mass loss through winds and jets 


Underestimated 


25% uncertainty on distance 


Error of a factor 2 


Inclination effects 


Underestimated 


Fitting procedure 


Error of a factor 2 


Radiative efficiency overestimated 


Underestimated 


Recurrence time overestimated 


Underestimated 



neutron stars), overestimating At by 50% will underestimate 
M ex t by a factor 1.5. 

A last point should be stressed regarding the recur- 
rence time. Many examples of transient sources show inter- 
outburst times varying from 1 to more than 10 years for 
the same binary system which can lead to significant errors 
if one uses an averaged value of At. Moreover, for several 
sources, the duration of an outburst is comparable or even 
greater than the time spent in quiescence (e.g. GX 339-4 
since 2002). Since the total mass accreted during an out- 
burst should correspond to the mass accumulated in the disc 
during quiescence plus the mass transferred to it during the 
outburst, neglecting the duration of the outburst in estimat- 
ing At can overestimate M ex t by a factor of 2. To minimise 
this, we calculated M ex t for each individual outburst by di- 
viding the total energy liberated during the outburst by the 
time elapsed between the end of the previous outburst and 
the end of the studied one. We then averaged the values of 
Me X t obtained for each individual outburst. 

Table [5] summarises the impact of the main sources of 
uncertainty on the estimates of the mass transfer rate. It 
appears from this table that a symmetric error of a factor of 
3 (quadratic sum of two factors of 2) should be taken into 
account but that, overall, Mext is likely underestimated. 



4 RESULTS AND DISCUSSION 

Figures Q] and plot the mass transfer rate against the or- 
bital period for the NS and BH sources respectively, together 
with the predicted separation between transient and persis- 
tent systems in the irradiated and non-irradiated cases. To 
plot the stability criteria we assumed a = 0.1 and C = 10 -3 . 
For the black hole systems, we assumed a primary mass Mi 
varying between 3M© and 15M© and for neutron stars we 
assumed Mi — 1.4Mq. Finally, for both types of systems we 
considered a mass ratio q varying between 0.1 and 1. Using 
the above parameter space, the stability criteria can be writ- 
ten Merit = fc-Ph r with, for the non irradiated case, b — 1.76, 
and k = (2.6 ± 0.9) x 10 16 for both the BH and the NS. In 
the irradiated case, b = 1.59, k = (3.9 ± 1.6) x 10 15 for the 
BH and k = (2.9±0.9) x 10 15 for the NS. The stability limits 
are represented by the shaded grey areas on Fig. [1] and [2] 

4.1 Neutron stars 

From a global standpoint, the neutron star data are in 
very good agreement with the theoretical predictions of 
the irradiated DIM. Indeed, we note that the stability line 
draws a fairly accurate limit between persistent and tran- 
sient sources. Assuming that the main uncertainty in deriv- 
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Cyg X-2 

SMC X-1 
LMC X-4 
Cen X-3 



Orbital period (hr) 



Figure 1. Mass transfer rate as a function of the orbital period for XRBs with neutron stars. The transient and persistent LMXBs have 
been indicated with filled and open symbols respectively, while the crosses indicate the high-mass, persistent systems. The shaded grey 
areas indicated 'DIM irr' and 'DIM non irr' represent the separation between persistent (above) and transient systems (below) according 
to the disc instability model when, respectively, irradiation is taken into account and when it is neglected. The horizontal dashed line 
indicates the Eddington accretion rate for a 1.4 Mq neutron star. The upper limits on the mass transfer rate are due to lower limits on 
the recurrence time except for 1A 1744-361 for which it arises from a upper limit on the distance to the source. The three left closed 
arrows do not indicate actual upper limits on the orbital period of SMC X-1, LMC X-4 and Cen X-3. They emphasise that the radius 
of any accretion disk in these three high-mass XRBs is likely to be smaller than the one derived from the orbital period since they likely 
transfer mass by a (possibly focused) stellar wind instead of (or in addition to) fully developed Roche lobe overflow. In the legend, the 
solid horizontal line separates transient and persistent systems. The dashed horizontal line separates low-mass and high-mass persistent 
systems. We can notice from this plot that the irradiated DIM seems to reproduce the transient-persistent division very well. 



ing the stability criterion lies in the irradiation constant C, 
we can use the narrow separation between the two groups of 
sources to roughly constrain its value. As we see on Fig. [T] 
for the value C = 10 -3 we initially assumed, the persistent 
sources 4U 2129+12 and Her X-1 lie already on the sta- 
bility line. If we increase the irradiation constant, the line 
shifts toward lower mass transfer rate and for C = 1CP 2 
it reaches the transient systems 4U 1608-52, GRS 1747-312 
and XTE J1751-305. Consequently, if the physical parame- 
ters assumed to derive the stability criterion are correct as 
well as our estimates of the mass transfer rates, the irradia- 
tion constant for NS XRBs should range between 10 -2 and 
10" 3 . 

4.1,1 Cyg X-2 

The only binary in our NS sample which does not agree with 
the DIM prediction is the persistent Z-source Cyg X-2 which 
should be transient according to its position in the diagram. 
However, the irradiated DIM predicts that the critical mass 
transfer rate corresponding to the large orbital period of Cyg 
X-2 should be an order of magnitude above the Eddington 
limit for a 1.4M0 neutron star. Since we derived M ex t from 
the (supposedly Eddington limited) X-ray luminosity, we 



might significantly underestimate the actual mass transfer 
at the outer edge of the disc. If the latter is truly an order 
of magnitude above MEdd, either the X-ray emitting flow 
becomes radiatively inefficient or outflows carry away 90% 
of the mass transferred from the donor star or the DIM 
has to be modified to properly describe super-Eddington 
accreting systems (for instance, to consider an Eddington 
limited irradiating flux). 

Another explanation would be that Cyg X-2 is in fact 
a transient source undergoing a very long outburst in a 
way similar to the well known black hole transient GRS 
1915+105. Long orbital period systems imply large accre- 
tion discs which can potentially accumulate enough mass to 
fuel a decades long outburst. To test this idea, we estimated 
the disc mass of Cyg X-2 befo re the onset of a hy pothetical 
outburst using equation (8) of lTruss fc Pond (2006) and the 
expression of th e critical den s ity of the cold neutral state 
(E~ it ) given in lLasota et al.l l|2008l . appendix A). We ob- 
tainec0 a disc mass Mdi ac = 2.5 x W 28 g. Assuming that the 

2 Note that, following iTruss fc Done] j2006h . we also integrated 
the surface density over a portion of the disc comprised between 
0.1-Rdisc an d ^f-disc but we rather assumed an average surface 
density of S ( ^. it (0.2i?, c ii sc ) instead of £~ ;t (0.1-f?disc) assumed by 
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matter is accreted onto the compact object at the Edding- 
ton rate, such a disc should be able to sustain a > 80 years 
long outburst. Given that Cyg X-2 was discovered 46 years 
ago, it could indeed be a transient source. 

4-1.2 Remarks 

It is also worth noting that the Z-sources of our sample 
(Sco X-l, Cyg X-2, LMC X-2 and GX 349+2) lie on or very 
close to the Eddington limit, in agreement with the known 
characteristic of this class of neutron stars XRBs. We note 
as well that the high-mass systems SMC X-l and LMC X- 
4 are accreting at a super-Eddington rate, which is also a 
known possible character istic of accretion onto X-ray pulsars 
(|Basko fc Sunvaevlll976| y 

4.2 Black holes 

The result of the test for the black hole sources is also in good 
agreement with the theoretical expectations. However, the 
low number of persistent sources, the uncertainties associ- 
ated to their nature and the numerous upper or lower limits 
due to unconstrained recurrence time or distance, weaken 
the significance of the test. Beyond this general assessment, 
we discuss below some individual cases that do not clearly 
fit into this global picture. 

4.2.1 IE 1740.7-2942 and GRS 1758-258 

These two persistent sources should be transient, according 
to their position in the diagram. As already discussed in this 
paper, their classification as persistent sources is unclear and 
given their large orbital period, decades-long outbursts are 
definitely not excluded. On the other hand, their orbital 
periods have been derived from modulations obser ved in 
their long-term RXTE lightcurves (|Smith et alj|2002l ) which 
might be less reliable than an optical determination. Con- 
sequently, we cannot exclude shorter orbital periods which 
could shift their position in the "persistent zone" of the dia- 
gram. For the above mentioned reasons, we won't consider 
these sources as contradictory to the DIM predictions. 

4-2.2 The size of the disc in Cyg X-l 

We note on Fig. [2] that Cyg X-l lies below the stability line 
of the irradiated DIM. However, as we mentioned previously, 
the size of its accretion disc might be smaller than the one 
derived through its orbital period. If the persistent/transient 
division is correctly described by the irradiated DIM, we can 
estimate an upper limit on the size of the disc of Cyg X-l, 
so that it is a persistent source given its mass transfer rate. 
Equating Eq. [2] with our estimated M ext for Cyg X-l and 
using a black hole mass of 14.8Mq we obtain an outer disc 
radius of 3.6 x 10 11 cm. This implies that the accretion disc 
of Cyg X-l should be at least 3 times smaller than the size 
derived assuming Roche lobe overflow. 

iTruss fc Done! 11200611 . We th ought this va lue was slightly more 
consistent with the results of lDubus et al.l (l200lh 



4.2.3 Is Swift J1753. 5-0127 a new persistent source? 

Swift J1753. 5-0127 was discovered by the Bur st Alert Tele- 
scope (BAT) on board Swift on 2005 May 30 jPalmer et all 
l2005h and, at the time of writing, is still undergoing the 
same outburst. No transition to the canonical soft state has 
been reported so far and the source always remained in hard 
states (Soleri et al. submitted). Given its particularly long 
outburst, its short orbital period (3.2h) and its position in 
the diagram of Fig. [5] (i.e. in the persistent zone) we investi- 
gated the possibility that Swift J1753. 5-0127 is a "new per- 
sistent system". We estimated an upper limit to the mass 
of the disc at the onset of the outburst by assuming that 
the surface density at all radii is equal to the critical sur- 
face density required to trigger an outburst- Given that the 
disc mass increases with radius, we assumed a 15 Mq black 
hole and obtained Af disc = 2.2 x 10 25 g. The 0.1-200 keV 
lightcurve integrated over time then gives us the total mass 
accreted from the start of the outburst until 2012 January 
1: Maccr = 6.6 x W 25 g. To estimate M acC r we u sed the lower 
limit to the distance to the source (7.2 kpc; IZurita et al.l 
2008) and a radiative efficiency n — 0.1. If the lower limit 
to the distance is correct, M accr should be considered as 
a lower limit to the total mass accreted since the fluence 
in the 0.1-200 keV range leads to a lower limit on the mass 
accreted. Consequently, even assuming that the entire (max- 
imum) disc mass was consumed in the outburst, the lower 
limit of the mass accreted so far is 3 times higher and the 
outburst is still ongoing. If we assume that the difference be- 
tween the mass accreted and the disc mass has been trans- 
ferred by the donor star along the outburst, we derive a mass 
transfer rate of M ex t ~ 2 x 10 17 gs _1 which is an order of 
magnitude above the critical mass accretion rate. Although 
very simple and crude, this calculation could indicate that 
Swift J1753. 5-0127 is a newly born persistent source. How- 
ever, considering the extremely low probability of observing 
such an event, it is likely that the mass transfer rate is in 
fact enhanced during the outburst by, e.g., irradiation of the 
donor star. If our calculation is realistic, the enhanced M oxt 
is greater than the critical accretion rate, which means that 
the outburst will last until the enhancement of the mass 
transfer weakens. 



4.3 A measure of the source activity - the 
transient ness parameter 

From a theoretical point of view, for a given orbital period, 
we would expect transient systems located closer to the sta- 
bility line to be more active, i.e. to undergo outbursts more 
regularly. Indeed, for a given size of the disc, the time needed 
to "refill" the disc after an outburst and reach the critical 
density necessary to enter in a new outburst will be lower 
if the mass transfer rate from the companion is higher. To 
test this idea we can define the "transientness" parameter 
of a source as the ratio between the estimated mass transfer 
rate and the critical accretion rate calculated for the orbital 



3 This assumption lead to an upper limit on the mass of the disc 
since the only requirement to trigger an outburst according to the 
DIM is that the surface density is above the critical value at a 
single radius only. 
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Figure 2. Mass transfer rate as a function of the orbital period for XRBs with black holes. The transient and persistent sources have been 
indicated with filled and open symbols respectively. The shaded grey areas indicated 'DIM irr' and 'DIM non irr' represent the separation 
between persistent (above) and transient systems (below) according to the disc instability model when, respectively, irradiation is taken 
into account and when it is neglected. The horizontal dashed line indicates the Eddington accretion rate for a 10 Mq black hole. All 
the upper limits on the mass transfer rate are due to lower limits on the recurrence time. The upper limits on the mass transfer rate 
of 4U 1957+115 and GS 1354-64 result from lower limits on the distance to the sources. The three left closed arrows do not indicate 
actual upper limits on the orbital period of Cyg X-1, LMC X-1 and LMC X-3. They emphasise that the radius of any accretion disk 
in these three high-mass XRBs is likely to be smaller than the one derived from the orbital period since they likely transfer mass by a 
(possibly focused) stellar wind instead of fully developed Roche lobe overflow. In the legend, the solid horizontal line separates transient 
and persistent systems. The dashed horizontal line stresses that the persistent nature of IE 1740.7-2942 and GRS 1758-258 is unclear 
(see Section 13) . 



period of the source. Table and 0] give the values of the 
transientness for our selected sources. 

The transientness of different sources can then be com- 
pared with their average recurrence time. From the DIM 
point of view, we expect the average recurrence time to de- 
crease with the transientness. However, the DIM, in its stan- 
dard form, predicts a regular recurrence time and identical 
outbursts, which is not observed in the majority of cases. In 
particular, some weak outbursts might be only echoes (i.e. 
XTE J1550-564, see below) of a "real" outburst as defined 
by the DIM. To estimate the average recurrence time of a 
given system, we therefore averaged the inter-outburst times 
weighted by the fluence of their subsequent outburst. This 
should minimise the bias introduce by short recurrence time 
leading to weak outbursts. 

Note that the black hole system XTE J1550-564 under- 
went a major outburst in 1998 (year of its discovery) fol- 
lowed by 4 weak outbursts, one every year, until 2003. Since 
then, the source remained in quiescence. Estimating the av- 
erage recurrence time is thus relatively complex in this case. 
In order to take into account the inter-outburst time pre- 
ceding the major outburst in 1998, we calculated the time 
necessary to accumulate the mass that was accreted during 
the 1998 outburst using our estimate of the mass transfer 



rate. We obtained ~ 28 years, which led to a weighted aver- 
age recurrence time of 20 years when including the 4 others 
inter-outburst times. 

Fig. OandUpresent the recurrence time as a function of 
the transientness for the neutron stars and the black holes 
respectively. We have only selected sources for which the re- 
currence time is constrained by at least two outbursts. The 
error bars are presented for indicative purpose and corre- 
spond to uncertainties of a factor 3 on the mass transfer 
rate and of 20% on the recurrence time. Given the large 
uncertainties, we cannot draw firm and quantitative con- 
clusions on the relation between the transientness and the 
recurrence time, but we note clear trends in these two fig- 
ures which are consistent with the DIM predictions: systems 
with higher transientness tend to undergo outbursts more 
frequently. This result strengthens again the DIM as the 
correct model to describe the global long-term activity of 
X-ray binaries. 



5 SUMMARY AND CONCLUSION 

We have revisited and refined a simple yet fundamental test 
of the disc instability model for X-ray binaries, first carried 
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Figure 3. Average recurrence time as a function of the transient- 
ness parameter (M C xt/M cr jt; see section|4]3jl for the neutron stars 
sample. 



points toward a mass transfer rate higher than the critical 
accretion rate. 

Finally, we introduce the transientness parameter as a 
measure of the activity of a source and show a clear trend of 
the average recurrence time to decrease with transientness 
in agreement with the DIM expectations. 

We therefore conclude that, despite some difficulties in 
reproducing the complex details of the light curves, the DIM 
succeeds in explaining the global behaviour of X-ray binaries 
averaged over a long enough period of time. 
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out by Ivan Paradiisl l|l996f ). Our work included a sample of 
52 sources and confirmed that the DIM is able to explain 
the dichotomy between transient and persistent sources as 
long as irradiation is taken into account. 

In the case of neutron stars, the data are in precise 
agreement with the DIM predictions which allowed us to 
constrain the irradiation constant between 10 -2 and 10 -3 . 
Although less constrained by the data due to the low num- 
ber of persistent sources and the numerous upper and lower 
limits, the black hole sources appear also in good agreement 
with the stability criterion predicted the DIM. 

We then discussed some individual cases which do not 
clearly fit into this general conclusion, and mainly showed 
that the classification of a source as a persistent system is 
arbitrary if the accretion disc is large enough to produce 
outburst longer than our historical record of its activity. We 
also emphasised the peculiar case of the black hole candidate 
Swift J1753. 5-0127 which might have already accreted more 
mass that the total maximum mass of its accretion disc. This 
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